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Abstract: We prove resolvent estimates for semiclassical operators such as — h 2 A + V(x) in scattering situations. 
Provided the set of trapped classical trajectories supports a chaotic flow and is sufficiently filamentary, the analytic 
continuation of the resolvent is bounded by h~ M in a strip whose width is determined by a certain topological pressure 
associated with the classical flow. This polynomial estimate has applications to local smoothing in Schrodinger 
propagation and to energy decay of solutions to wave equations. 
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1 Statement of Results 



In this short note we prove a resolvent estimate in the pole free strip for operators whose classical Hamiltonian 
flows are hyperbolic on the sets of trapped trajectories (trapped sets), and the latter are assumed to be sufficiently 



<**i ^ filamentary - see (|1.4[) for the precise condition. The proof is based on the arguments of (|2ll ) and we refer to §3 



of that paper for the preliminary material and assumptions on the operator. 

The polynomial estimate on the resolvent in the pole free strip below the real axis 1|1.5|) provides a direct 
proof of the estimate on the real axis (|1.6p . and that estimate is only logarithmically weaker than the similar 
bound in the non-trapping case (that is, the case where all classical trajectories escape to infinity). Through an 
argument going back to Kato, and more recently to Burq, that estimate is crucial for obtaining local smoothing 
. and Strichartz estimates for the Schrodinger equation. These in turn arc important in the investigation of 
nonlinear waves in non-homogeneous trapping media. Also, as has been known since the work of Lax-Phillips, 
the estimate in the complex domain is useful for obtaining exponential decay of solutions to wave equations (see 
the paragraph following (|1.6|) for some references to recent literature). 
. An example of an operator to which our methods apply is given by the semiclassical Schrodinger operator 

O ; Pu(x) = P(h)u(x) = -h 2 — V d Xj (V$9 ij d Xi u(x)) + V(x) , xeW 1 , (1.1) 

OS . V9 ^ 

o 

G(x) d = (g 1 -' (x))i,j is a symmetric positive definite matrix representing a (possibly nontrivial) metric on W 1 , 

g d = 1/ detG(a;), and V(x) is a potential function. We assume that the geometry and the potential are "trivial" 
outside a bounded region: 

g lj (x)=Sij, V(x) = -1, when|x|>i?. 

This operator is hence associated with a short-range scattering situation. We refer to (|2ll . §3.2) for the complete 
set of assumptions which allow long range perturbations, at the expense of some analyticity assumptions We note 
that for V = — 1, P(h)u — is the Hclmholtz equation for a Laplace-Beltrami operator, with h = 1/A, playing 
the role of wavelength. 

Such operators have a purely continuous spectrum near the origin, and their truncated resolvent 
x{P{h) — z)~ x x (x € C£°(1R")) can be meromorphically continued from Imz > to Imz < 0, with poles of 
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finite multiplicity called resonances. In the semiclassical limit the distribution of resonances depends on 

the properties of the classical flow generated by the Hamiltonian 

n 

that is the flow (x,£) i— » exptH p (x,£) associated with the Hamiltonian vector field 

n 
k=l 



(when V = — 1 the Hamiltonian flow corresponds to the geodesic flow on S*M. n .) More precisely, the properties 
of the resolvent x(P(h) — Z )~ 1 X near z = are influenced by the nature of flow on the energy shell {p(x, £) = 0}. 
A lot of attention has been given to nontrapping flows, that is flows for which the trapped set 

K d = {(x,£) : p{x,i) = 0, exptH p (x,€) ~h oo, f -> ±00} (1.2) 

is empty. In that case, for 5 > small enough and any C > 0, the resolvent is pole free in a strip [— 6, 5] — i[0, Ch], 
and satisfies the bound (jl7l : 

- «)- 1 xlU»-£» - OQT 1 ), z G [-5, <5] - i[0, C7l] . 

On the opposite, there exist cases of "strong trapping" for which the trapped set has a positive volume; resonances 
can then be exponentially close to the real axis, and the norm of the resolvent be of order e c ^ h for z 6 [—5, 5] 

In this note we are considering an intermediate situation, namely the case where the trapped set (|1.2p is a 
(locally maximal) hyperbolic set. This means that K is a compact, flow-invariant set with no fixed point, such 
that at any point p 6 K the tangent space splits into the neutral (M.H p (p)), stable (EZ~), and unstable 
directions: 

T p p-\0) = RH p (p) 8 E~ ffi E+ . 

This decomposition is preserved through the flow. The (un)stable directions are characterized by the following 
properties: 

3A>0, \\dexptH p (p)v\\ < C e~ A|t| ||u|| , Vn € Ej, ±t > 0. 

Such trapped sets are easy to construct. The simplest case consists in a single unstable periodic orbit, but we 
will rather consider the more general case where K is a fractal set supporting a chaotic flow; such a set contains 
countably many periodic orbits, which are dense on the set of nonwandering points NW(K) C K (fl5h ) . 

Our results will depend on the "thickness" of the trapped set, defined formulated in terms of a certain 
dynamical object, the topological pressure. We refer to (21. §3.3) and texts on dynamical systems (jl5l : 12!^ for 
the general definition of the pressure, recalling only a definition valid in the present case. Let / 6 C°(K). Then 
the pressure of / with respect to the Hamiltonian flow on K is given by 

V{f) = r Hm^log Y, ex P / \^PtH p rf( Pl )dt, (1.3) 

where the sum runs over all periodic orbits 7 of periods T 7 < T, and p 7 is a point on the orbit 7. The function 
/ we will be using is a multiple of the (infinitesimal) unstable Jacobian of the flow on K: 

<P+(p) d = -17 det(dexptH p \ E+)\ t=0 , p € K . 
We can now formulate our main result: 

Theorem. Suppose that P(h) satisfies (jl.ip or the more general assumptions of (00, $3.2). Suppose also that 
the Hamiltonian flow is hyperbolic on the trapped set K , and that the topological pressure 

V(—ip+/2) < 0, (p+ the unstable Jacobian. (1-4) 



Then for any x € C£°(K™) and e > 0, there exist 5(e) > and h(e) > smc/i that the cut-off resolvent 
X(P(h) — z)~ l X, Imz > ; continues analytically to the strip 

Q e (h) d = {z : Imz > h(V(-tp+/2) + e) , \Rez\<S(e)}, < h< h(e) . 

For z G f2 e (/i) H {Imz < 0}, this resolvent is polynomially bounded in h: 

\\x(P(h) - z^xW^, < C(e, X ) h~ 1+CE lmz/h log(l// l ) , 

dcf n (1.5) 

CE 2|P(-^+/2) + e/2|- 

For any s <E [0, 1], the pressure T'(—sip + ) measures relative strengths of the complexity of the flow on K 
(i.e. the number of periodic orbits), and the instability of the trajectories (through the Jacobian). For s = 0, 
V(0) only measures the complexity, it is the topological entropy of the flow, which is generally positive. On 
the opposite, V(— <p+) is negative, it represents the "classical decay rate" of the flow. The intermediate value 
T'(~(p+/2) can take either sign, depending on the "thickness" of K. In dimension n = 2 the condition (|1.4p is 
equivalent to the statement that the Hausdorff dimension of K C p -1 (0) is less than 2. Since the energy surface 
p _1 (0) has dimension 3 and the minimal dimension of a non-empty K is 1, the condition means that we are less 
than "half-way" and K is filamentary. Trapped sets with dimensions greater than 2 are referred to as bulky. 

The first part of the theorem is the main result of (f2lh . see Theorem 3 there. Here we use the techniques 
developed in that paper to prove (|1.5[) . For the Laplacian outside several convex obstacles on K n (satisfying 
a condition guaranteeing strict hyperbolicity of the flow) with Dirichlet or Neumann boundary condition, the 
theorem was proved by Ikawa (UJ), with the pressure being only implicit in the statement. For more recent 
developments in that setting see (|2Tl . ((iil ) . and (0). 

In particular, for z on the real axis the bound (| 1 . 5j) gives 

\\ X (P(h) - z^xh'^L' < C l0g{ V h) , z€ [S(e),S(e)}, < h< h(e) . (1.6) 

This result was already given in (|2ll . Theorem 5) with a less direct proof. It has been generalized to a larger 
class of manifolds in (|9|) and fl 1 .5(1 provides no new insight in that setting. 

One of the applications of (11. 6p in the case of the Laplacian is a local smoothing with a minimal loss (0) in 
the Schrodinger evolution (see |j) for the original application in the setting of obstacle scattering) : 

VT>0,Ve>0,aC = C(!7», f \\xe{-itA g )u\\ 2 H1/2 ^dt<C\\u\\ 2 L2 . 

Jo 

One can also deduce from (|1.6p a Strichartz estimate (0; HI) useful to prove existence of solutions for some related 
semilinear Schrodinger equations. 

In the case of the Laplacian (V = —1), the estimate in a strip (|1.5p has important consequences regarding 
the energy decay for the wave equation - see (0; S E3) and references given there. In odd dimension n > 3, 
it implies that the local energy of the waves decays exponentially in time. The same type of energy decay 
(also involving a pressure condition) has been recently obtained by Schenck in the setting of the damped wave 
equation on a compact manifold of negative curvature (j2~il ). 

To prove (|1.5p we use several methods and intermediate results from |2lh. Using estimates from (|2~ll §7), 
we show in f}3]how to obtain a good parametrix for the complex-scaled operator, which leads to an estimate for 
the resolvent. As was pointed out to us by Burq the construction of the parametrix for the outgoing resolvent 
was the, somewhat implicit, key step in the work of Ikaw a (|14ri on the resonance gap for several convex obstacle. 
That insight lead us to re-examine the consequences of (|2ll). 

We follow the notation of (1211) with precise references given as we go along. For the needed aspects of 
semiclassical microlocal analysis (21, §3) and the references to (fioh and (fill ) should be consulted. 



2 Review of the hyperbolic dispersion estimate 

The central "dynamical ingredient" of the proof is a certain dispersion estimate relative to a modification of 
P(h), which we will now describe. 

The first modification of P{h) comes from the method of complex scaling reviewed in (|2lL §3.4). For any 
fixed sufficiently large Rq > 0, it results in the operator Pg(h), with the following properties. To formulate them 
put 

n e d £ [-5,6] +i[-e/C,C] , = M 1 h\og{l/h). (2.1) 

Then 

Pg(h)-z : Hl(R n ) — ► L 2 (W l ) is a Fredholm operator for z £ tt e , (2.2) 
VxeC c °°(S(0,i?o)), X R(z,h)x = xRe(z,h)x. (2.3) 

Here and below we set the following notation for the resolvents: 

R.(z,h) d = {P.(h) - z)- 1 , lmz>0, 

and (|2.3p shows the meromorphic continuation of %i?(z,/i)x to Qg, guaranteed by the Fredholm property of 
Pe{h)-z. 

The operator Pg(h) is further modified by an exponential weight, G w = G w (x, hD), 
GeC™(T*W l ), suppGcp" 1 ((-2 < 5,2,5)), d a G — C(/ilog(l//i)) , 
where S > is a fixed small number. The modified operator is obtained by conjugation: 

P e>e (h) = e- eGW ' h P e {h)e eGW ' h , e = M 2 6, 6 = Mi/ilog(l//i) . (2.4) 

This operator has the same spectrum as Pe(h) and the following properties: 

if tpa £ S(T*W"-), supp^o Cp- 1 ((-3,5/2, 3J/2)), , 

if Imip%{x, hD) P e ,e{h) $f(s, hD) <Ch (2.5) 

(2.6) 

The main reason for introducing the weight G is to ensure the bound (|2.5p . The specific choice of G is 
explained in (j2ll . §6.1). In particular G vanishes in some neibhbourhood of the trapped set K, and the operator 
exp(eG w (x, hD)) is an ft,-pseudodifferential operator B w (x, hD), with symbol satisfying 

Beh- N S 5 (T*W>), B ] CsuppG =l + Os 5 (h°°). 

As a result, if the spatial cutoff x is supported away from 7rsuppG, calculus of semiclassical pseudodifferential 
operators ensures that 

xR(z,h)x = xReA^h)x + L 2^(h 00 )\\Re, e (z,h)\\. (2.7) 

From now on our objective will then be to estimate the norm ||i?g ie (z, h)\\ L 2^ L 2 . 

We consider a final modification of Pq^(K) near the zero energy surface. Let ipo € S(T*M. n ) be supported 
in p _1 ((— 3(5/2, 35/2) and equal to 1 in p~ 1 (—5,5). Define 

P , e (h) = ffl(x, hD) P e , t W(x, hD) , (2.8) 

and the associated propagator 

U(t) d = exp{-itP e , e {h)/h} . (2.9) 
The crucial ingredients in proving (|1.5[) are good upper bounds for the norms 

\\U(t)ip w (x,hD)\\ L 2^ L 2 , on time scales < t < M log(l//i) , 



where M > is fixed but large, and 

V>eS(T*R"), suppV> Cp- 1 ({-6/2, 5/2)), ip = 1 on p -1 ((-<y/4, 5/4)) . (2.10) 

From the bound (|2.5[) on the imaginary part of Pg^ e (h), we obviously get an exponential control on the propagator: 

\\U(t)\\ L 2^ L 2 < exp(Ci), i>0. (2.11) 

The reason to conjugate Pg with the weight G w was indeed to ensure this exponential bound. Together with the 
hyperbolic dispersion bound (|2. 13[) . this exponential bound would suffice to get a polynomial bound 0(h~ L ) in 
Q1.5|) . for some (unknown) L > 0. To obtain the explicit value, 

c e Im z 
+ _ ~ ' 

for the exponent, we need to improve (|2.f f [) into the following uniform bound: 

Lemma 2.1. Let tp satisfy the conditions (|2.10p . Then, there exist ho, C'o > smc/i t/iai, 

||i/(t)t/) u '(af ) ftZ>)|| I ,a_ > i,2 < C , < t < Mlog(l//i), /i</i . (2.12) 

Before proving this Lemma, we state the major consequence of our dynamical assumptions for the classical 
flow on K , namely its hyperbolicity and the "filamentary" nature of K (expressed through (If .4jl ). ft is a hyperbolic 
dispersion estimate which was explicitly written only in a model case (|2lL Proposition 9.1) (|2ll ). but can be easily 
drawn from (f2ll . Proposition 6.3), in the spirit of (21, §6.4) As above, we take tp as in (|2.10p . For any e > we 

set A = f —V((p+/2) + e/2. For any < h < h(e), we then have 

\\U(t)i> w (x,hD)\\ L ^ L2 < Ch- n / 2 exp(-\t)+0(h M *), 
uniformly in the time range < t < M\og{\/h) . 

The constant M is arbitrarily large, and M3 can be taken as large as we wish, provided we choose M\ in (|2.1[) 
large enough depending on M. If the pressure V((p+/2) is negative, one can take e small enough to ensure 
A > e/2 > 0. The above estimate is then sharper than (I2.12|) for times beyond the Ehrenfest time 

t E d ^c E \og(l/h), c E d ^^. (2.14) 

The large constant M will always be chosen (much) larger than ce- 

Proof of Lemma \2.12\ To motivate the proof we start with a heuristic argument for the bound (|2.12j) . As 
mentioned above, the exponential bound (|2.1ip is due to the fact that the imaginary part of Pg^(h) can take 
positive values of order 0(h) (see (|2.5jl ). However, the construction of the weight G shows that outside a bounded 
region of phase space of the form 

V pos =p-\(-25,25))nT{ Ri<lxl<R2} R\ 

the imaginary part of Pg^(h) is negative up to 0(h°°) errors. 

The radius i?i above is large enough, so that V pos lies at finite distance from the trapped set. As a result, 
any trajectory crossing the region Vp OS will only spend a bounded time in that region. For this reason, the 
propagator U(t) on a large time t > 1 will "accumulate" exponential growth during a uniformly bounded time 
only. 

We now provide a rigorous proof, using ideas and results from (|2ll . §6.3). The phase space T*R™ is split 
using a smooth partition of unity: 



1= nb > 7T b eC°°(T*R n ,[0,l}). 



b=0,l,2,oo 

These four functions have specific localization properties: 



• supp7r 6 C Sj) for b = 0, 1, 2 

• 7TOO is localized outside 35/4, 35/4)) 

• 7Ti is supported near K, in particular, its support does not intersect V^, s 

• 7T2 is supported away from K but inside {\x\ < R2 + 1} 

• 7To is supported near spatial infinity, that is on {|x| > R2 — 1} where the operator Pg. e (/i) is absorbing (the 
imaginary part of its symbol is negative) . 

Employing a positive (Wick) quantization scheme (see for instance flit), and for the semiclassical setting 
((23I §3.3)), Tib = OPh (""&)) we produce a quantum partition of unity 

Id= n &< HUfcll < 1 - 

b=0,l,2,oo 

The evolution U(t) is then split between time intervals of length to, where to > is large but independent of h. 
Using the partition of unity, we decompose the propagator at time t = Nto into 

U(Nt )ip w (x,hD) = ( U b ) N ip w (x,hD), where U b d = U(t ) IL b . 

fc=0,l,2,oo 

Expanding the power, we obtain a sum of terms U bN ■ ■ ■ U bl ip w ; to understand each such term semiclassically, 
we investigate whether there exist true classical trajectories following that "symbolic history", namely sitting 
in supp^ at time 0, in supp7rf, 2 at time to, etc. up to time Nto. 

Since the energy cutoffs ip and tToo have disjoint support, no classical trajectory can spend time in both 
supports. As a result, any sequence containing at least one index bi = 00 is irrelevant (meaning that the 



corresponding term is Ol2^ L 2 (h°°)) (|2ll . Lemma 6.5) 



Since any classical trajectory can travel in supp7r2 at most for a finite time < iVoto before escaping, ((21 
Lemma 6.6) shows that the relevant sequences 61 - - ♦ 6jv are of the form 

bi = l for N Q < i < N - N . 

They correspond to trajectories spending most of the time near K. One then has 

U(Nto) i> w {x, hD) = U(Noto) (Uif- 2 ^ U(N t ) ip w {x, hD) + Q L ^ L i{h M ^) , 

uniformly for any 2 No < N < Mlog(l//i), where M5 > is large if the previous M, Mi are. 
Finally, using the fact that the weight G vanishes on supp7ri, |2ll Lemma 6.3) shows that 

u 1 = u{to)n l = Uo{to)n 1 + o L ^^{h 00 ), 

where C/o(to) — cxp(—itoP{h)/h) is unitary. Hence, \\Ui\\ < 1 + 0(h°°), while ||[/(Afoto)|| is estimated using 

(gup . □ 

3 Resolvent estimates 

We can now prove the resolvent estimate (|1.5j) by constructing a parametrix for Pg i€ (h) — z, z G fl e (h) defined 
in the statement of the theorem. We will use the notation 

C = z/h 

to shorten some of the formulae. We want to find an approximate solution to 

(P g Jh)-z)u = f, feL 2 (R n ), zen e (h). 

First, the ellipticity away from the energy surface p _1 (0) shows that, for ip as in (|2.10p . there exists an operator, 
T =0(1) : i 2 (R") H%(R n ), such that 

(P B ,e(h) ~ *W = (1 - i> w {x, hD))f + R f , R - L ^ L *(h°°) . 



To treat the vicinity of p 1 (0) we put 

T 1 f=(i/h) dte i(t U{t)^ w (x,hD)f, t M = Mlog(l/h), 



Jo 

which satisfies 

(P e , e (h) - z) Tif = ip w (x, hD)f + Rif , Ri = -e^ U(t M ) ^ w (x, hD) . (3.1) 

The estimate (|2.13[) shows that, if A + Im£ > e/2, and for arbitrary M4 > 0, one can choose M and M3 large 
enough such that R\ = Ojji^ip. (h Mi ). We can estimate the norm of T\ by the triangle inequality, 

< h- 1 / e- Im(t \\U{£)il> w (x,hD)\\ L 2^ L 2dt, (3.2) 
Jo 

and then use the bounds (|2.12p for times < t < Ie and (|2. 13|) for times t& < t < tjvf. 

When ImC = 0, the above integral can be estimated by the integral over the interval t € [0, tg]: 

ImC = 0=^ HTxIU^^ < h- 1 (c t E + j S ) < Ch- 1 \ogh- 1 . 

In the case > lm( > —A + e/2, the dominant part of the integral comes from t = tg: 

> ImC > -A + e/2 WT^^^^ < CjT r e~ lniCtE = C e h- 1+CE ImC . 

We rewrite Q3.ip as 

W(x, hD)(P e<e {h) - z)^(x, hD)Tif = i> w (x, hD)f + R t f . 
From the inclusion V'olsuppi/j = 1, one can show (as in (|2lL Lemma 6.5)) that 

i/)%(x t hD)(Pe te (h)-z)il>%(x,hD)Ti = (P e , e (h) - z)Ti + R 2 , R2 = L ^ L .{h°°) , 

and also that 

llTill^ClllilU,. 

h 

Putting T = T + T 1 and R = R + Ri + R 2 , we obtain 

(Pe, e (h) - *)T = ld+R, R = L 2^ L2 (h m ) . 
This means that (Pg te (h) — z) can be inverted, with 

\\(PeAh) - zy'h^Hl - (1 + 0{h M *))\\T\\ L2 ^ Hl . 

h h 

The above estimates on the norms of To and T\ can be summarized by 

0>ImC>e + ^(-y+/2)=^ \\T\\ L ^ H 2 < C £ h- 1+CE ImC log/i" 1 . (3.3) 

Using (P?rp . this proves the bound (|1.5p . □ 

Remark. By using a sharper energy cutoff tph belonging to an exotic symbol class (see |27l . §4)) and supported 
in the energy layer p _1 ((— ft, 1- " 5 , h 1 ^ 6 )) (as in (1)), the bound (|2.13p is likely to be improved to 

\\U(t)^(x,hD)\\ L 2^ L 2 < Ch- {n - 1+S)/2 e* v {-\t) + 0{h Ah ) . (3.4) 

This bound becomes sharper than (|2.12p around the time t' E = c' E \og(l/h), where 

. def n — 1 + S 

C E = — jj— < CE. 



As a result, the bounds on the norm of the corresponding operator T[ are modified accordingly. At the same 
time, as shown in yj, Prop. 5.4), the ellipticity away from the energy surface provides an operator Tq satisfying 

(Pe, e (h) - z)U = (1 - ffi(x, hD)) + L 2^ L 2 (h°°) , 

and of norm ||T^|| i 2^ jL 2 = 0(h- 1+s ). The norm of T = T^+T[ is still dominated by that of T{, so that we 
eventually get 

\\ X (P(h) - z^xWl^hI < Ce h- 1+c '* lmz ' h log(l//i) , z e n e (h) n {Imz < 0}. 



Since it is not clear that even this bound is optimal, and that proving (|3.4|) would require some effort, we have 
limited ourselves to using the established bound (|2.13p . 

One advantage of the approach presented in this note (compared with the method of (|2ll . §9)) is that, to 
obtain the bound (|1.6p we did not have to use the complex interpolation arguments of Q) and |28h. 
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